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' Abstract. We present a stability result for a wide class doubly nonlinear equations, featu- 

i ring general maximal monotone operators, and (possibly) nonconvex and nonsmooth energy 

functionals. The limit analysis resides on the reformulation of the differential evolution as a 
i-O ■ scalar energy-conservation equation with the aid of the so-called Fitzpatrick theory for the re- 

presentation of monotone operators. In particular, our result applies to the vanishing viscosity 
MH ' approximation of rate-independent systems. 

\ 1. Introduction 

' This note is concerned with a convergence result for doubly nonlinear differential inclusions of 

, the type 

a n (u n (t)) +d£ t (u n (t)) 3 in 1* for a.a. t e (0, T). (1.1) 

Here, (a n ) is a sequence of maximal monotone (and possibly multivalued) operators a n : X =4 X*, 
(X, || ■ ||) is a (separable) reflexive Banach space, and £ : [0, T]xl^ (— oo, oo] is a (proper) time- 
dependent energy functional. We will prove that for a n — > a in the graph sense any limit point u 
\ of the sequence (w n )neN is & solution to 

a(u(t)) +dE t (u(t)) 3 in X* for a.a. t e (0, T) 



OO 



o 



Throughout the paper, we write £t(it) in place of £(i,w). We will understand the multivalued 
operator <9£ : (0,T)xI^t X* to be a suitable notion of subdiffercntial for the possibly nonsmooth 
and nonconvex map u i-> £t(u), namely the so-called the Frechet subdifferential, defined at (t, u) € 
dom(£) by 

£ G <9£((m) if and only if > £*(u) + (£,v — u) + o(\\v — u\\) as v — > u. (1.2) 

Observe that, as soon as the mapping it i— > £t(u) is convex, the Frechet subdifferential 9£ t (u) 
coincides with the subdifferential of it h- > £t(u) in the sense of convex analysis. 

Doubly nonlinear equations as in (jl.ljl arise in a variety of different applications, ranging from 
Thcrmomechanics, to phase change, to magnetism. As such, they have attracted a substantial deal 
of attention in recent years. Correspondingly, the related literature is quite rich. Being beyond 
our scope to attempt here a comprehensive review, we limit ourselves to recording the seminal 
observations by Moreau [56l [57] and Germain [37], as well as the early existence results by 
Arai 0, Senba [71] . Colli & Visintin [26], and Colli [25] . The reader can find a selection of 
more recent results in [SJ [21 O [JJ [33l Ell EU [65l ES]. Without going into details, let us mention 
that, over the last decade, the convexity requirement on the map u H ► £ t (u) in the pioneering 
papers [9[ [7TJ [26l [25] has been progressively weakened: in particular, in [51] a quite broad class 
of nonsmooth and nonconvex energy functionals has been considered. Nonetheless, in all of the 
aforementioned contributions the operator a is assumed to fulfill some coercivity property, namely 
to have at least linear growth at infinity. We will refer to this case as viscous. 
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The case of O-homogeneous operators a has been recently investigated as well, for it connects 
with the modeling of so-called rate-independent systems. We shall hence refer to this situation as 
rate-independent. Some references in this direction are to be found in the papers [28, 29, 36] 1431 

BS EH [51 M HI]. 

Additionally, relation has been considered in connection with the study of the long-time 
behavior of solutions in [H [2U [SHI EH US] , and their variational characterization in [5] [5] (72] [75] . 

The focus of this paper is on the study of the stability of the doubly nonlinear flows (jl.ip . 
Namely, we investigate the convergence of solutions u n to equations under the assumption 

a n — > a in the graph sense in X x X* , (1-3) 

viz. for all (£, £*) £ graph(a) there exists (£ n , £*) 6 graph(a„) such that — > £ in X and £* — ► £* 
in X* as n — > oo. The main result of this paper, Theorem 14.51 states that cluster points u of the 
curves (it„) are in fact solutions to the limiting equation 

a (ii{t)) + d£ t (u(t)) 3 in X* for a. a. t £ (0, T). (1.4) 

We have to mention that stability results for the doubly nonlinear flows are already 

available in the literature. For viscous graphs a n , a first convergence theorem in the case of 
convex energies has been obtained by Aizicovici & Yan [1] (see also [72]), whereas stability 
results for doubly nonlinear equations with nonconvex energies have been proved in [51] . This 
issue has been recently reconsidered by Visintin [801 HU H2] , who has remarkably extended the 
reach of the theory to treat subdifferential inclusions of the type 

£(«(*)) + 7(u(t)) 3 inX* for a.a. t £ (0, T), 

with /3, 7 : X =t X* maximal monotone operators, /3 cyclically monotone and 7 noncyclic mono- 
tone, by resorting to the so-called Fitzpatrick theory [34] . 

Let us briefly recall that an operator a : X =4 X* is cyclically monotone if a is the generalized 
gradient of some potential. Namely, if a = dip for some proper, convex, and lower semicontinuous 
function tp : X — > (—00,00], where the symbol d here denotes the subdifferential in the sense of 
convex analysis. In the cyclic-monotone case a = dip, it is well known that the relation y £ dip(x) 
can be equivalently reformulated as (y,x) = tJj(x) + tl)*{y) 1 where ip* is the Legendre-Fenchcl 
conjugate of ip and (•, •) is the duality pairing between X* and X. The use of this variational fact 
for the aim of variationally reformulating evolution equations dates back to Brezis-Ekeland 
[151 1X7] and Nayroles [59j 160] . Among the many contributions stemming from this idea, the 
reader is especially referred to the existence proofs by Auchmuty [11] and Roubicek [67], and 
to the recent monograph by Ghoussoub 39J on self-dual variational principles (see also the 
references in |72)V 

The Fitzpatrick theory allows us to extend this variational view to subdifferential inclusions 
of the type (|1.4|) . with a possibly noncyclic monotone, by introducing representative functions 
f a :Xx X* — > (—00, 00] for the operator a. These are convex functions f a with the property 

M{x,y) £ X x X*, (y,x) < f a (x,y) and 

y £ a(x) iff (y,x) = f a (x,y). 

The reader is referred to Section [2] below for a selection of relevant results within this theory. 
In particular, in 80] these tools are used in order to reformulate variationally relations (|1.1[) 
for noncyclic monotone operators. This reformulation opens the way to devise a suitable T- 
convergence analysis toward structural stability of the flows. 

As for the rate-independent case, one shall mention the stability results for hysteresis operators 
from the classical monographs [HI [42] [78] (see also [74]). Another stability result in the rate- 
independent setting is in [72]. Moreover, we record Visintin [81, 82 , which exploits the Fitzpatrick 
idea in the rate-independent context, but by taking perturbations in d£ n (again, in a possibly 
noncyclic monotone frame). 

Finally, the approximation of rate-independent flows by viscous flows (in the cyclic-monotone 
case) has been recently attracted a great deal of attention. This is especially critical as viscous and 
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rate-independent evolutions usually call for different analytical treatments. The vanishing viscosity 
approach to abstract rate- independent systems has been in particular developed in [32] [3S] [5U] . 
More specifically, in the latter two papers it has been shown that the vanishing viscosity limit leads 
to the notion of BV solution to a rate-independent system. In the recent [52], still within the 
cyclic-monotone frame, the p n — > 1 limit, where p n is the homogeneity of the potential %p n of a n , 
has been addressed, and it has been proved that BV solutions arise in the limit. A stability result 
with respect to variational convergence for the latter solution concept has also been obtained. 

Our result. The focus here is that of obtaining a stability result for the differential inclusions 
([Lip by allowing for maximal generality on the perturbations a n and on the functional £. In 
particular, we shall neither assume super-linear equi-coercivity in a n , nor cyclic monotonicity. As 
for the energy £, we do not require neither smoothness nor convexity with respect to u, but still 
ask for lower semicontinuity and some coercivity, see Assumption 13.91 below. 

This generality sets our result aside from the available contributions on this topic. In particular, 
our analysis also encompasses the passage from viscous to rate-independent doubly nonlinear 
evolution. Indeed, we are able to treat here the p n — ¥ 1 case for noncyclic monotone operators a n 
(in this setting, p n is the coercivity exponent for a n ). In our general context, we prove that the 
so-called local solutions [49J [50] to a rate-independent system arise in the p n — > 1 limit. 

The basic idea for handling the noncyclic monotone case, is to resort to a variational reformu- 
lation of the flows (|l.ip which is well suited for discussing limits. By letting f an represent the 
monotone operator a n in the sense of (|I.5I) and assuming the validity of a suitable chain rule 
for the energy £, relation (|l.ip is proved to be equivalent (see Proposition I3.I2[ ) to an energy 
conservation identity, namely 



The strategy is then to prove that, by passing to the liminf in (|I.6p . the structure of the relation 
is preserved. In particular, we provide sufficient conditions under which the liminf of the integral 
of the representative functions f an is a representative function of the limit graph a. Care here 
is given to developing such a lower semicontinuity argument for functions which are only BV in 
time. This allows us to directly include the case of rate-independent flows. 

Let us once more emphasize that we can encompass in our analysis a broad class of time- 
dependent energies £ : [0,T]xI-> (— oo, +oo], (possibly) nonsmooth and nonconvex with respect 
to the state variable it, but still satisfying a suitable set of coercivity and regularity type conditions 
mutuated from [51] . Note that in such a general setting existence of absolutely continuous solutions 
u n : [0,T] — > X to (|l.ip is presently not known. A possibility in order to overcome this would 
be to strengthen our assumptions on the energy functional £, for instance assuming it to be a 
suitable perturbation of a (A-)convex functional, as in [55] [551 \5T\ . We however refrain from this, 
for the sake of keeping maximal generality for the convergence result. In this respect, our result 
should be regarded purely as a stability analysis, with focus on the convergence properties of the 
operators (a n ). A stability result with respect to suitable convergence of the energy functionals 
could also be obtained, closely following the lines of [52 Thm. 4.8]. Again, we have chosen not 
to detail this, in order to highlight the usage of the Fitzpatrick theory to deal with the noncyclic 
operators (a n ). This very generality will allow us, for instance, to address in the upcoming Section 
14.31 the quasistatic limits of a class of dynamical problems, which can be in fact reformulated as 
doubly nonlinear equations of the form (jl.ip . 

Structure of the paper. Section[2]contains some background material on Fitzpatrick theory, and 
on the notions of variational convergence for functionals and operators which will be relevant for 
the subsequent analysis. In Section [3. 11 some further preliminaries of measure theory and convex 
analysis are provided, whereas in Section 13.21 the basic assumptions on the energy functional £ 
are stated in detail, and suitable reformulations of (II. ip are discussed. In Section 2] we state our 
main stability result Theorem 14.51 and thoroughly comment it. We also give two corollaries (i.e. 




energy at t v v 

dissipated energy on [0, t] 



initial energy v v 

work of ext. actions 




(1.6) 
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Theorem 14.81 and Proposition 14. 9|) in two particular cases: specifically, Prop. 14.91 deals with the 
p n — > 1 vanishing- viscosity limit. We conclude Sec. U by discussing classes of energy functionals 
to which our results apply (cf. Sec. I4.2[) . and developing applications to rate- independent limits of 
Hamiltonian systems (in Sec. I4.3[) . The proof of Theorem 14.51 is developed throughout Section [3 
also exploiting some results from Young measure theory which are contained in Appendix El 

2. Fitzpatrick theory 

Within this section, we shall systematically use the notation 

tt&O := (CO for all G A x X* 

for the duality pairing between the reflexive space X and X*, and identify possibly multivalued 
operators a : X =$ X* with the corresponding graphs a C X x X* without changing notation. 
We recall that a : X =4 X* is monotone if 

<t for all C G a(0, Co G «($>) 

(where (•, •) denotes the duality pairing between X* and X), and that it is maximal monotone, if 
it is maximal for set inclusion within the class of monotone operators. 

We shall provide here a minimal aside on the Fitzpatrick theory, essentially mutuated from [79 . 
The reader is however referred to I2U 1231 13H1 IM1 HS1 HS1 Ell E3] for additional material and 
a collection of related results and applications to PDEs. 

Representative functionals. We denote by 7(A) the set of functionals <p : X x X* —> (—oo, oo] 
such that 

<p is convex, lower semicontinuous, and y>(£,£*) > 7]"(£,£*) V(£, £*) G X x X*. 
We associate with </? G ^(A) the set a C X x X* given by 

(e,r)Ga ^ <p(z,n (2.i) 

Whenever (|2.1I) holds we say that ip represents a, that cp is representative, and that a is repre- 
sentable. A representable operator can be represented by different representative functionals (cf. 
Example 12 .41 below) . On the contrary, each representative functional represents only one operator. 

Example 2.1. In the cyclically monotone case of a — dip for some (proper) convex and lower 
semicontinuous potential ip : X — > (—00,00], a representative functional for a is given by the 
bipotential (according to the terminology of [20j ) 

^,O = V>(0 + V>*(O- (2-2) 

We have the following strict set inclusions 

{maximal monotone operators} 5 {representable operators} S {monotone operators}. 

Namely, representable operators are intermediate between monotone and maximal monotone. One 
may wonder how to translate maximality at the level of representative functionals. The following 
result provides a useful criterium for the representability of a maximal monotone operator. 

Proposition 2.2 (Representative of a maximal monotone operator [76]). A functional <p G ^(X) 
represents a maximal monotone operator iff ip* G 5"(X*). In this case, if ip represents a then ip* 
represents cT 1 = {(£,£*) : (£*,£) G a}. 

The Fitzpatrick and the Penot functions. Given a C X x X* with a ^ I we define the 
Fitzpatrick function (associated with a) /„ : X x J* ^ (-00, 00] by 

= sup{7r(£ ,n- t(£o-£,£) : (&>,«)€ a} for all (CD G^ x A* 
and the Penot function (associated with a) p a : X x A* — > (—00, 00] by 

p a := (tt + /„)**. (2.4) 
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Both f a and p a represent a. Moreover f a and p a are respectively the minimal and maximal 
element (with respect to pointwise ordering) of the functional interval 

3(a) = {ip e 5(X) : ip represents a}. (2.5) 

In particular, in the cyclically monotone case of a = dip, there holds 

ut,n<m+r(n<p a (Ln fa ex xx: (2.6) 

Let us also point out that, in view of the definitions (|2.3j) and (|2.4|) of f a and p a the following 
formulae hold 

/afen=sup{ 7 r^o,n+^^o)-Pa- I (^^o) ■ (6,S)6Jfxr} V e XxX* , (2.7) 
Pa(^r)=sup{7r(^o,n+^^o)-/a-(^»Co) : (£o, £o) e XxX*} ¥(£,£*) eXxX*, (2.8) 

namely f a (p a , resp.) is the convex conjugate of the Penot (Fitzpatrick, resp.) function of the 
inverse operator a -1 . 

Finally for later use we observe that 

Oea(0) >0 for all (£,£*) e A x A*. (2.9) 

Self-dual representatives. Recall that a function tp : X x X* — > (—00, 00] is called self-dual iff 

= ^(r,0 for all (CD G A x X*. 

The interval J(a) from (|2.5|) includes self-dual representative functions [631 Thm. 3.3]. In the 
cyclically monotone case of a = dip, an example in this direction is given by the bipotential 
(I2.2[) . Out of the cyclically monotone realm, an example of a self-dual representative in the case 
a = 7 + dip with 7 skew adjoint is (£, £*) 1-^ -0(C) + V>*(-7£ + C) El]. 

In the general case, the indirect proof of the existence of self-dual representative functions is 
due to Penot [6TJE2] an d Svaiter [76], whereas direct constructions have been firstly provided 
by Penot & Zalinescu [53] under some restriction on a. An explicit self-dual representative 
function in the general maximal monotone case has been recently obtained by Bauschke & Wang 
[15] and reads 

(cn^J i n t v + c, r + o + - c, r - o + iicii 2 + ncii2}. (2.10) 

Note that neither the Fitzpatrick function f a nor the Penot function p a are self-dual in general. 

Let us now recast the characterization of maximal monotonicity of Proposition 12.21 in the fol- 
lowing. 

Proposition 2.3 (Self-dual representatives = maximality). An operator a : X =4 X* is maximal 
monotone iff it is represented by a self-dual functional (p. 

Proof. By [TS] , if a is maximal monotone, then it admits the self-dual representative (|2.10l) . 
As for the converse implication, note that by self-duality of <p and it we get 

o = ¥>&o >^,r) = 7r*(r,o- 

Thus, p* e J(X*) and Proposition 12.21 applies. □ 

Self-dual representatives vs. Fitzpatrick and Penot functions. Let a C X x X* be a 

cyclically monotone operator with a = dip for some convex and lower semicontinuous potential 
ip : X — >• (—00,00]. As already observed, a self-dual representative of a is the sum of ip and its 
convex conjugate. However in general the Fitzpatrick functional f a may differ from ip + ip*, as 
shown by the following. 

Example 2.4. Consider X = M = X* and set a = identity, namely a = dip with ip(£) = \£, 2 - 
The Fitzpatrick function of dip is / s ^(£, £*) = £ 2 /4 + (£*) 2 /4 + £ ■ TA which is not self-dual. 
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Fitzpatrick and Penot functions in the case of 1-positively homogeneous potentials. 

Our next result reveals that, when a = dip and ip : X — > (— oo, +00] is positively homogeneous of 
degree 1, then also the Fitzpatrick functional f a coincides with the bipotential (|2.2p . 

Proposition 2.5. Let ip : X — > (—00, 00] be convex, lower semicontinuous, and positively homo- 
geneous of degree 1, viz. ip(\£) = A^>(C) for all C £ X and A > 0. 

Then, the Fitzpatrick function of the subdifferential of ip coincides with the sum of ip and its 
convex conjugate, i.e. 

fo4tO=m + r(€) forall(Z,Z*)£XxX*. (2.11) 

Before developing the proof we recall that, for all ip convex, lower semicontinuous, and 1-homogeneous 
there exists a closed convex set G K C X such that ip coincides with the Minkowski functional 
of K, viz. 

V»(0 = M K (£) = inf la > : - £ k\ . (2.12a) 



(2.12b) 



Furthermore, 

ip*{C) = Ik* (C*) where K* c X* is the polar set of K, i.e. 
K* = {CEK* : <!V£eA-}, 
so that 

V({,0^ M K (Z) + I K .(?) (2.12c) 
For later convenience, we also recall that 

ip(£)= sup tt(C,C*) foraU^eX. (2.12d) 



Proof of Proposition \2.5[ It follows from the definition of the Fitzpatrick function (|2.3j) and from 
(|2.12c|) that 

fwfoC) = sup{7r(&,e*) - Mk(Co) +tt(^^) - Ik* (Co*) : (Co, Co) G <*} 

V v ' " v ' 

<m k (s) 

Hence, we obtain that /»0 (£,£*) < ^it(C) + ^if*(C*) for all (C,C*) £ X x X*. 

For the opposite inequality assume first that £* ^ K* . Then there exists Co £ -K" such that 
^(CcC*) > 1- Choose an arbitrary Co £ <9Mr-(Co) = 9Mk(ACo) for any positive A > 0. Then, 
taking into account that M^(ACo) + ^*(Co) = -^(CojCo) — ^> we S et 

t(ACo,C*) - A/^(ACo) + tt(C,Co) - Ik4Co) 

> A(tt(Co,C*) - 1) +^(C,Co) ^ +°o as A — > 00. 

Therefore /a^.(C,C*) > ^k*(C*)- On the other hand, taking into account that dM K (0) = K*, we 
deduce that 

iW(C,C*) > sup{7r(C,C *) : Co* e iT} = A%(C) . 

Eventually, we get that /a^(C,C*) > m k{£) + Ik*(€*) for all (C,C*) £ X x X*, which concludes 
the proof. □ 

Corollary 2.6. Let ip : X — > (—00, 00] be convex, lower semicontinuous and positively homoge- 
neous of degree 1 . Then 

hAt,n = m+r(n = Pd<P'(^o foraii(s,oexxx*. (2.13) 

Proof. To prove (|2.13p we observe that 

/av = ^ + V>* = (i> + ip*T = fdi, = po4>- 

where the first identity is due to Proposition ^. 51 the second one to the fact that ip + ip* is self-dual, 
and the last one to (12.81). □ 
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2.1. Approximation of maximal monotone operators. In the next lines, the symbol X may 
stand for the space X, for X*, or for X x X*. Let f n , f : X — > (—00, 00] be convex, proper, and 
l.s.c. functionals and let a n , a C X x X* be maximal monotone operators. We introduce the 
notation 

r-liminf f„ (x) := min{liminf f n (x„), x n — x in X}, 

n— J-oo n— >oo 

T-limsup f n (x) := min{limsupf„(x„), x„ — > x in X}. 

n— >oo n— >oo 

These correspond to the classical r-liminf and r-limsup constructions (cf. e.g. [17]), with respect 
to the weak and the strong topology of X, respectively. We will use the following convergence 
notions, for which the reader is referred to |10) : 

f„ — » f <^ r-limsup f n < f < r-liminf f„ 

n— >oo n—too 

a n a V((,(*)ea 3(( n ,(*)6a B : £n £*■ 

where the symbol stands for Mosco convergence in X and — is usually referred to as graph 
convergence. In particular, Mosco convergence corresponds to T-convergence with respect to both 
the strong and the weak topology of X and can be made more explicit by 

Vx„ x, f(x) < liminf f„(x n ) 

n— too 

Vx S X, 3x„ -> x : f„(x„) ->• f(x) 



f„ > f <S> 



In the case of cyclically monotone operators, graph convergence is known to be equivalent to 
the Mosco convergence of the respective potentials (up to some normalization condition), viz. we 
have the following result. 

Theorem 2.7 ([T01 Thm. 3.66, p. 373]). Let <fi n , <f> : X — > (—00, +00] be proper, convex and 
lower semicontinuous functionals. The following are equivalent: 

i) d(j) n — — > d<f> in X x X* and there exist (£n,£n) £ d^n such that — ^ £ in X, £* — > £* in 

X*, 0(0, and(t,?)ed4>; 

ii) <j) n — > 4> in X . 

The importance of graph convergence is revealed by the following identification lemma, basically 
consisting of the approximation version of [16| Prop. 2.5, p. 27]. 

Lemma 2.8. Let a n a, (£ n ,C) e a n, £n ^ Q ^ liminf^oo tt(£„, £*) < 7r(£,D- 
Tfterc (£,£*) G a. 

In order to prove this, quite classical, approximation lemma, in the particular case of cyclically 
monotone operators, what is actually needed is just the implication i) =4> ii) in Theorem l2.7l above. 

Indeed, let d<j> n d(j>, (£„,£*) G 5 0„, with £n 6 C ^ C, and liminf T ,_ ! , O0 tt(C„,C) < t(S,D 
be given. We readily have that 

< 0(0 + 0*(f) - Trfef) < liminf (</>„(£„) + - = 



n— 7-oc 



where we have used the fact that 0* — — >• 0* iff </>„ — >• 10, Thm. 3.18, p. 295], i.e. the 
bicontinuity of the Legendre-Fenchel transformation with respect to the topology induced by the 
Mosco convergence. Therefore we conclude that (£,£*) E a. 

Remark 2.9. As indeed one just needs 

d(f> n d(p =S> (f> < r-liminf <j> n in V and cf)* < r-liminf (j>* in V* 

n— >oo n— too 

in order to check for Lemma 12.81 one may wonder if asking directly the two r-liminf conditions 
above would weaken the convergence requirements on the functionals. This is however not the 
case. Indeed, under some very general condition of equi-properness type, we have that the two 
separate r-liminf conditions are indeed equivalent to (f> n — >■ (f> \72\ Lemma 4.1] and hence entail 
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Our next aim is that of extending the above arguments to the case of noncyclically maximal 
monotone operators. In particular, we present an extension of Theorem 12.71 in terms of represen- 
tative functions, and in particular of the Fitzpatrick and of the Penot functionals. 

Theorem 2.10. Let (a n ), a be maximal monotone operators a n , a : X =4 X* . The following are 
equivalent: 

i) a n — —> a in X x X* , 

H) fa< r -liminf f a in X x X* , 

n— >oo 

Hi) r-lirnsup n _ s . 00 p an < p a in X x X* . 

Exactly as above, the proof of Lemma [2.81 follows just from implication i) ii). We shall 
however give a full equivalence proof for the sake of completeness and comparison with Theorem 
12.71 In particular, note that Condition ii) above is weaker than f an f a . That is to say that 
the former Theorem 1 2 . 71 does not follow directly from Theorem 12. 101 



Proof. Claim 1: i) => ii). Fix (£o,£o) e a an d ^ (£on,£on) G a n be such that £o« ~> Co m X 
and £q„ — > £q in X*. Moreover, let £„ — »■ £ in X and £* — £* in X*. We have that 

liminf / Q „(£n,Cn) > liminf (7r(£ n,O ~ ^(Con - £n,£oJ) 

n— >oo n— >oo 

= ^o,n-^o-^,e )- 

In particular, by passing to the supremum with respect to (£o>£o) € a, we conclude that f a < 
T-liminf n^oo f an . 

Claim 2: i) => Hi). Observe that a n — — > a implies a" 1 — — > a -1 , hence Claim 1 yields f a -i < 
T-liminf / -i. By convex conjugation and taking into account (12.81) and [101 Thm. III. 3. 7, p. 271], 

n— >oo ™ 

we then have 



P« = (/a-0* > (r-hminf / -i = r-limsupp« n . 

Claim 3: m) => i). Fix (£,£*) G a and let (£ n ,££) £ a « fulfill £„ ->■ £ in X, £* -> in X*, and 

lim sup p Q „ < p Q (£,£*) = 

(such sequences exist as r-limsup^^^ p Qjl < p a ). In particular, we have that 

for some sequence e n — » 0. By exploiting the extension of the Br0nsted-Rockafellar approximation 
Lemma from [44J Thm. 3.4], we have that there exist (£m£n) € ot n such that for all n G N 

Then, a classical diagonal-extraction argument yields — > £ in X and £* — > in X*. 
Claim 4: u) =>• i). Again by convex conjugation, and (|2.8I) . we deduce from that 

Pa-i > r-limsupp Q -i. 

n— too 

Therefore, in view of Claim 3 we have that aC 1 — — » a -1 , whence a„ — — > a. □ 



3. Setup and preliminary results 

Before stating our working assumptions in Sec. l3.2[ in the upcoming Sec. l3.ll we recall all the basic 
definitions, and tools of measure theory and convex analysis, which we will use in the following. 
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3.1. Preliminaries of measure theory, BV functions, and convex analysis. We start with 
the notion of measure with values in a Banach space X, which later on will either coincide with 
the reflexive space X, or with R. 

Definition 3.1 (Vector measure). Let (r2,£) be a measurable space. A function /i : E — >• X is 
called a (Banach-space valued) vector measure, if 

V (Ai) ieS , A, e E with (i^j^ At n Aj = 0) it holds fi I (J A t J = £ M (At) . (3.1) 

YieN / i=l 

Here the convergence of the series on the right-hand side has to be understood in terms of the 
norm of X. 

Definition 3.2 (Variation of a measure). Let (f2, S) be a measurable space and fj, : £ — > X a 
vector measure. Then the variation of //, denoted by \\fx\\ : S — )• [0, oo], is given by 

M(A) := sup j 53 11^)11 = (-Ai)ieN C £, |J ^ = A , Vi ^ i : ^ n A,- = 1 (3.2) 

li=l i<EN ) 

for all A £ £. If ||/z||(f2) < oo then we say ^ is of bounded variation. 
Indeed, ||ju|| itself is a (positive) measure on (0,S), see [301 Prop. 9, p. 3]. 

Definition 3.3 (Absolute continuity and singularity of measures). Let (0,E) be a measurable 
space, (i:E->Xa vector measure, and v : E — > [0, oo] be a (real- valued, positive) measure. We 
say that /i is absolutely continuous w.r.t. to v, and write /i <C if 



Vies 



=0 =^> /Lt(A) = . (3.3) 



Moreover, we say that two real-valued, positive measures fi and v are singular, and write [i _L v, 
if there exist Si, B 2 G E with Si U B 2 = ^ and Si n B 2 = such that 

VAeE: = (i(An B±) and = v (A n B 2 ) (3.4) 

We recall the following generalization of the Lebesgue decomposition theorem, see e.g. [30l Thm. 
9, p. 31]. 

Theorem 3.4 (Lebesgue decomposition theorem). Let (0,E) be a measure space, a be a Banach 
space-valued measure of bounded variation and X a real valued, positive measure. Then there exist 
two unique vector measures cr aC; a s - m on (f2, E), which are of bounded variation, such that 

Uracil < A, ||cr sin || _L A and a = er ac + cr sin • (3.5) 

BV functions. We fix here some definitions and notation concerning BF-functions on [0, T] with 
values in a Banach space X, referring e.g. to |58j for a comprehensive introduction to the topic. 
We denote by BV([0, T];X) the space of the measurable, pointwise defined at every time t € [0, T], 
functions v : [0,T] — > X such that their pointwise total variation on [0,T] is finite, i.e. 

Var(u; [0,T]) = sup j ^ ||v(t TO ) - u(t ro _i)|| : = t < *i < ... < t M -i < t M = t| < oo . 

More in general, given a convex, lower semicontinuous, 1-positively homogeneous functional ip : 
X — > [0, +cxd), we denote by Var^, the induced total variation, i.e. 

Vai>(w; [0, T]) = sup j 53 ^(«(*m) ~ u(*m-i)) : = t < h < . . . < t M -i < t M = t\ < oo . 

(3.6) 

It is well known that the distributional derivative dv of a curve v £ BV(0,T; X) is a vector 
measure in M(0,T;X), where 

M(0, T; X) = { Radon vector measures fj, : (0, T) — > X with bounded variation }, 

which we will endow with the weak* -topology. 
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Notation 3.5. Let u e BV([0,T];X). Applying Thm. I3H with the choices a = du and A = £ 
(where £ denotes the one-dimensional Lebesgue measure on [0, T]), we find that there exist vector 
measures (du) ac , (du) s i n 6 M(0, T;X) such that 

||(dtt) ac |] < £, ||(du) sin || 1 £ and du = ( du)^ + ( du) sia . (3.7) 

Thanks to the Radon-Nikodym property of the reflexive space X, the Radon-Nikodym derivatives 

Uac(t) == ~77T— , Uri»(t) := j^i 51 " exist for a.a. t e (0, T). (3.8) 
d£ IKdujginll 

For later use, we remark that for any convex, lower semicontinuous, 1-positively homogeneous 
if> : X — > [0, +oo ) there holds 

Var^(u;[0,r])= / *P(u ac (t))dt+ [ ^{u sin {t))\\{du) sin \\{t) . (3.9) 



The recession function. Finally, we recall the concept of recession function (see [351 Chap. 4]). 
Note that the following definitions and results, which are stated in [35J for convex functions on 
R m , in fact extend to an infinite-dimensional setting, as it can be easily checked. 

Definition 3.6 (Recession function). Let X be a vector space and g : X — > (— oo, oo] be a convex 
functional. Its recession function g°° is defined as 

g°°[z) := sup {g(y + z) - g(y) : y e D(g)} . (3.10) 

Trivially adapting the argument from [35J Thm. 4.70, p. 290], it can be shown that that g°° is 
positively homogeneous of degree 1 and convex. Moreover, if g is lower semicontinuous, so is g°° . 
Furthermore, there holds 

oo, \ r g(y + tz)-g(y) g(y + tz)-g(y) . . 

g (z) = hm = sup for every y S -D(g) . (3.11) 

In what follows, we will denote by the recession function of the Fitzpatrick function f a , viz. 

/~&0=™P + +x*)~fa(x,x*) : (x,x*) e D(f a )} . (3.12) 
We now prove a useful representation formula for f™, cf. [35, Prop. 4.77, p. 294]. 
Lemma 3.7. There holds 

/~(e,n = su P {(r,eo) + <€s,o : (M)^^-)} neneixr. (3.13) 

Proof. Following the proof of [351 Prop. 4.77], from (|3.1ip and (|2.7p we infer 

/a Is i ? ) — bu P 7 
t>0 t 

t>0 t V 

+ (£a> x ) + ( x *,£o) - Pa-^Co^o) ~ fa(x,X*fj 

> <C,£o) + (M 

+ \ «4M + (s*,6>) - - /a(a:,a:*)) Vfo.ft) e A x A*, i > 0. 

In view of (I2.7|) . we thus conclude that 

£°(£,0>sup{<r,£o> + <£o.0 : (fo,Qe%-i)}. 

The converse inequality may be proved arguing along the very same lines, cf. also the proof of [35i 
Prop. 4.77]. □ 

As a direct consequence of Lemma 13.71 we have the following representation formula for the 
recession function of f a , in the case a is the subdifferential of a 1-positively homogeneous potential. 
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Corollary 3.8. Let ip : A — > K be convex, lower semicontinuous and positively homogeneous of 
degree 1, and let K* C X* be the associated polar set, cf. (|2.12p . Then, 

f^,C)=sM(C ,O + {C,^} : (to,e o )eX xK*}. (3.14) 

Proof. Formula (|3.14j) follows from (|3.13[) . taking into account that 

Pa- 1 = Pdi>* = t/) + ip* = t/> + I K , , 
and that D(tp) = X by assumption. □ 

3.2. Basic assumptions. In what follows, we will suppose that 

X is a reflexive Banach space (3.15) 

and that 

a : X =$ X* is a maximal monotone operator with 6 a(0). ([3]ao) 

As for the energy functional £, along the lines of |51] we require the following coercivity and 
regularity type conditions. Recall that d£ denotes the Frechet subdifferential of the map u n- 

£t(u), cf. 02). 



Assumption 3.9 (Assumptions on the energy). We assume that the pair (£, d£) has the following 
properties: 

Lower semicontinuity: The domain of £ is of the form D(£) = [0,T] x D for some Del, 
and d£ : [0,T] xDzj X* . Furthermore, we ask that 

u i V £ t (u) is l.s.c. for all t e [0,T], 3C > : V(f,it) e[0,T]xfl: £ t (u) > C and m . 

graph(5£) is a Borel set of [0, T] x X x X* . ^ o) 

Coercivity: Set 5(u) := sup te [ t] £t( M ) f or every u 6 D. We require that 

u i-> 9 (it) has compact sublevels. (|3j £ i ) 

Time- differentiability: For any u 6 D the map t i— > £j(it) is differentiable with derivative 
dt£t(u) and it holds 

BdXl-.VueD: \d t £ t {u)\ < Ci£ t (u) ©£3) 

Weak closedness: For all t G [0,T] awe? /or sequences (u n ) n< £iq C A, £ n e d£t(u n ), E n = 
£ t (u n ) and p n = d t £ t {u n ) with 

u n — > u in A, — ^ ^ m A*, p„ — >• p, anc? E n E in M. 

it holds 

(t,u) <E D(dE), £ £dE t {u), p<d t £-t{u) and E = £ t (u) . (0l£ 3 ) 

Remark 3.10. In fact, up to a translation, we may always suppose that the constant involved in 
(|31£n[) is strictly positive. As in [BTj . combining (1 3 1 £ ? [) with the Gronwall Lemma we observe that 



3C>0 V(t,ii)e[0,r]xfl 5(u)<C inf £ t (u). (3.16) 

te[o,T] 



Later on, Assumption 13.91 will be complemented by a suitable version of the chain rule for £, cf. 
Assumption 14.41 below. As already mentioned, in order to investigate the stability properties of 
the doubly nonlinear equation 

a{u{t)) + d£ t {u{t)) 3 in A* for a.a. t € (0, T) , (3.17) 

under graph convergence of a, it is essential to resort to the Fitzpatrick function f a associated with 
a. In the following lines, we will therefore shed light on how p,17[) can be in fact reformulated in 
terms of an energy identity (cf. ()3.20[) below) featuring f a . At first, we will confine the discussion 
to the case of absolutely continuous solutions u to (|3.17|) . 
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Reformulations of (|3.17l) in the absolutely continuous case. Preliminarily, let us precisely 
define what we understand by an absolutely continuous solution to (|3.17[) . 



Definition 3.11 (Absolutely continuous solution). In the framework of ()3.I5|) . pictnj) . and (|3l£nj) . 

we say that a curve u G M /1,1 (0, T; A) is a solution to (|3.I7p . if there exists £ 6 i 1 (0, T; A*) with 

£(*) G (-a(u(t))) n d£ t (u(t)) for a.a. t G (0, T) . (3.18) 

In what follows, with a slight abuse of notation we will sometimes say that (u, £) is a solution to 
(|3~T7|) . meaning that (|3TT51) holds. 

In Proposition 13 . 1 2l we reformulate (13.181) by means of an energy identity involving the Fitz- 
patrick function f a . In the proof, a key role is played by the chain-rule condition (|3.19j) below on 
the energy £, whereas note that not all of the conditions collected in Assumption 13.91 are needed. 

Proposition 3.12 (Variational reformulation). In the framework of (|3.15[) . let a : X =4 X* fulfill 
(|3la:n[) and suppose that £ : [0,T] x X — > (-co, +oo] complies with (|31£n|) . (|31£i [) . (|31 £-^>|) . and the 
following chain rule: for every u G W l,1 {0,T; X) and £ G L 1 (0,T; X*) such that 

sup £ t (u(t)) < oo, £(t) G dE t {u{t)) for a.a. t G (0,T), / f a {u{t), -£(*)) di < oo, 
te[o,T] Jo 

(observe that, thanks to (|31 the first of the conditions above guarantees \dt&t(u(t))\ dt < oo 
as we/?), toere ZioWs 

</ie map 1 1— > £ t (tt(i)) is absolutely continuous and 

d (3.19) 
— £ t («(t)) - <£(t), u(t)> + 9e£t(«(t)) V 0.0. te (0, T). 

Then, the following implications hold: 

(1) if{u,£) G W /1,1 (0, T; A) x L^OjTjA*) /uZ/iZfe i/ie energy identity 

£t(u(t))+ / / a (u(«), -£(*)) ds=£ (u(0))+ / d t £ s (u(s))ds /or aZI i G (0, T], (3.20) 

i/ien (it, £) is a solution to (I3.17P in £/te sense of Def. \3.11\ 

(2) every solution (it,£) to (|3.17j) (m i/te sense of Def. [3.11]) fulfilling 

sup £*(«(*)) < oo, / | (£(t),u(t)} \dt < oo, (3.21) 
te[o,T] Jo 

complies in addition with the energy identity (|3.20[) . 

Observe that, for every solution (u, £) to (|3.17l) . since — £ G a(w) a.e. in (0, T) and G a(0), we 
have (— £, u) > a.e. in (0, T). Hence the second of (13.21[) in fact reduces to J (— £, it) di < oo. 

Proof. Let (it, £) fulfill (|3.20[) . Taking into account that f a (?i, — £) > a.e. in (0,T) thanks to 
(|2.9|) . and exploiting (|3I£9[) we gather 

£*(«(*)) < £o(t»(0)) + Ci / £ s (n(s))ds for alU G (0,T], (3.22) 

Jo 

whence sup i£ r T i £t(u(t)) < oo. Therefore, a fortiori (|3.20p yields f a (u, — £) G i 1 (0,T). Hence 
the pair (it,£) fulfills the conditions for the chain rule (|3 . 19[) . which yields for all t G (0,T] 



(3.23) 



Using that f a represents a, it is immediate to deduce from the above inequality that — G 
a (m(£)) for almost all t G (0, T), thus (u, £) is a solution to f|3.17[) in the sense of Def. 13.111 



f en (u(s), -£(*)) ds < £ (u(0)) - £ t («(t)) + / a t £ s (u(s)) ds 
< / (-£(«), u(t)) ds. 
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Conversely, let (u, f) G W 1 - 1 ^, T; X) x L x (0, T; X*) be a solution to (j3~T71) (in the sense of Def. 
13. lip fulfilling in addition (|3.21j) . Then, since f a (u, — £) = (— the chain rule f|3.10[) applies, 
yielding, for all £ G [0, T], the energy identity 

r f a (u(s), -as)) ds = ( (-£(«), u(t)) ds 



= £ («(0)) - £ t (u(t)) + f &&M*)\ ds • 
Jo 

a 

Remark 3.13. A few comments on Proposition 13.121 are in order. 

(1) It is not difficult to check that in Proposition 13.121 the Fitzpatrick function f a could be 
replaced by any representative functional for a. 

(2) Observe that, in the chain of inequalities (|3.23[) leading to the proof of part (1) of Propo- 
sition [37T21 it is in principle necessary for (13.201) and for the chain rule (|3. 19[) to hold as 
inequalities, only. The proof of part (2) requires (|3.19j) to hold as an equality, instead. 

4. Main results 

Before stating Thm. l4~5l let us precise our hypothesis on the sequence (a n ) of maximal monotone 
operators. 

Assumption 4.1. Let a n : X X* fulfill (|3lanP for all n G N and 

3ci,C2,C3>0, p > 1, q>l VneN V(x,y)Ga„ : 

(y,z)> Cl ||:rr + c 2 |M|*- C 3. ^ 

Furthermore, there exists a : X =4 X* fulfilling (|3lan[) such that a n a. 

Remark 4.2. Combining piai [) with the graph convergence of (a„) to a, it is immediate to conclude 
(y,x)> Cl \\x\\ p + c 2 \\y\\l~c 3 for all (x, y) G a , (4.1) 

The following example guarantees that our analysis encompasses the p n — > 1 vanishing- viscosity 
limit. 

Example 4.3. Let (p n ) C [1, +oo) fulfill p n 4- 1 as n —> oo, and let us set 

^ n (a:) = — ||x|| p ™, a n = 0Vn 

Clearly, (-0„) Mosco-converges to ^(a;) = ||x||, hence (a n ) converges in the sense of graphs to 
a = dip. Observe that tp^iv) ~ S^IMI*" with q n — p n /{Pn — 1) G [2, oo] for all n G N, and that 

(y.z) = + = = llyllf. for all (x,y) G a„ . 

p n q n 

Therefore, Assumption 14 . 1 1 is satisfied. 

The main result of this section addresses the passage to the limit as n — ► oo in the doubly 
nonlinear equations 

a n {u{t)) + dZ t (u(t)) 3 in X* for a.a. t G (0, T) . (4.2) 

In particular, we will assume to be given a sequence (u n ) of absolutely continuous solutions to (|4.2[) 
and we will show that, if the sequence (a n ) complies with Assumption 14. ll up to a subsequence 
(ti„) converges to a curve u fulfilling a suitable generalized formulation of (I3.17[) . 

Observe that, (|31o;i [) in principle only allows for a bound of the type ||un||ii(o,T;X) < C. That 
it why, we can only expect a BV([Q, T]; X)-regularity for the limiting curve u, and (|3 . 1 T[) has to 
be weakly formulated accordingly. This will be done through an energy inequality akin to ()3.20|) . 
cf. (14. 5p below. Therein, suitable replacements of the "time-derivative" of u are suitably handled 
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in terms of the Fitzpatrick function f a and of its recession function (cf. Definition 13. 6[) , and 
of the absolutely continuous and singular parts of the Radon derivative du of u. Having in mind 
the role of the chain rule (|3.19[) relating (|3 . 1 T[) and the energy identity (|3.20|) , it is to be expected 
that a suitable BV version of (|3 . 19[) will play a relevant role. We state it in the following: 

Assumption 4.4. Let u G BV([0,T];X) and £ G L X {Q,T; A*) fulfill 

sup £t(u(t)) < oo, £(*) G d£ t (u(t)) for a.a. t G (0,T) , f /„(«(*),-£(*)) di < oo, 
te[o,r] Jo 

and suppose that the map t i— > £ t (u(t)) is almost everywhere equal on (0, T) to a function E G 
SV([0,T]). Furthermore let du and d£ denote the Radon derivatives of u and E. 

Then, for almost all Lebesgue points t of the absolutely continuous parts w ac and E ac of du 
and dE there holds 

E ac {to) > <£(to), Uac(*o)> + 9t£ t0 (u(t )) /or £(i ) G 5£ to (u(i )) • ©£4) 

Observe that, since X has the Radon-Nikodym property, the set of Lebesgue points of u ac and 
-Eac has full Lebesgue measure in (0, T). 

As it will be clear from the proof of Thm. 14.51 below. Assumption 14.41 does not only provide a 
motivation for the energy inequality f|4.5[) , but it also has a key role in the proof of the passage to 
the limit as n — > 00 in (14.21). 



Theorem 4.5. Assume ()3.15j) . Lei a„, a : A" =4 A* fulfill Assumption 14. ll and suppose that 
£ : [0, T] x A — > (— 00, +00] complies with Assumptions 1331 and l4~4l Lei ms consider a sequence 
(uq ) C D of initial data such that 

u'^u inX, £oK)) -> £oOo), (4-3) 

andfei (u„,f„) C PF 1 ' 1 (0, T; A) x L 1 (0, T; A*) &e solutions to (g^J) in the sense of Definition\MH 
fulfilling the initial conditions u„(0) = Mq. Suppose that, in addition, for all n G N rte functions 
(« n ,^„) comply with the energy identity (|3.20[) . 

Then, there exist functions u G BV([0,T];A) and £ G L 9 (0,T;A*) fimi/i q > 1 /rom (|3lcn [) ) 
satisfying u(0) = Uq, £(i) G <9£t(u(i)) /or almost allt G (0,T), and swc/i i/iai up to a (not relabeled) 
subsequence 

u n (t)^u(t)Vte[0,T], du n = (u n ) ac -£\ [0 , T ] A d7ieM(0,T;A), (4.4) 
and (it, £) satisfies the energy inequality 



(4.5) 



as we/Z as 

e(t)€(-a(«ac(*)))nfl£ t («(t)) for a.a.t G (0, T) . (4.6) 

Furthermore, there exists E G SV([0, T]) suc/i £/ia£ 

E(t) = E t {u{t)) for a.a.t € (0, T), E(t) > E t (u(t)) for all t G [0,T], (4.7) 

and we aave i/ie pointwise energy identity 

Kc(*)+ /a («ac(*), "£(*)) =ft£t(«(*)) /or a.a. i G (0, T). (4.8) 

Remark 4.6. In view of Proposition l3~T2l a sufficient condition for functions (u n ,£„) solving (|4.2I) 
to comply with the energy identity (|3.20p . is that they fulfill 

sup E t (u n (t)) < co and (-(„>„) ei'fO.T). 
te(o,T) 

This, provided that the absolutely continuous version ()3.19[) of the chain rule holds. 



£*(«(*))+/ /a(Wac(s),-£(s))ds+ / /~ (u sin ( S ), 0) [| ( d«) sin || (s) 

Jo Jo 
< £oO(0)) + / 9 t £ s (w(s))ds /or all t G [0,T], 
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In Section l4"T2l we will discuss some sufficient conditions on £ for both the chain rule p. 191) and 
its BV- version of Assumption 14.41 to hold. 



4.1. Further results. We conclude this section with some results which shed light on the inter- 
pretation of the energy identities (|4.5I) and (|4.8I) satisfied by the pair (u, £). More precisely: 

- Proposition ^. 71 focuses on the case in which we have the additional information that u is 
absolutely continuous. For instance, this is granted whenever u occurs as limiting curve 
of a sequence (u n ) C M /1,1 (0, T : X) of solutions to the differential inclusions (|4.2I) . driven 
by operators (a n ) which fulfill a stronger version of condition (|3lai [) . cf. Thm. FOl ahead. 

- In Proposition 14.91 we address the special case in which a = dip, with ip : X — > [0, +oo) 
a convex, lower semicontinuous, and 1-homogeneous dissipation potential. We show that 
in this case any u G BV([0, T]; X) complying with the energy inequality (|4.5[) is a local 
solution (cf. [1H1[SD]) to the rate-independent system (X, £.,ip). 

The absolutely continuous case. Under a slightly stronger version of the chain rule of As- 
sumption 14. 4| Proposition 14.71 shows that, if in addition we have that the curve u is absolutely 
continuous on (0, T), then (u s i n (t), 0) = for |[ ( dii) S i n |[-a.a. t G (0, T), and (|4.5|) holds on every 
sub-interval [s,t] C [0, T]. Furthermore, the pair (u, £) solves (|3.17j) in the sense of Definition 13. Ill 
cf. (|4~TT1) below. 

Proposition 4.7. In the framework of ()3.15j) . let a : X =4 X* fulfill (|3lan|) . and £ : [0,T] x 
X — ► (—00, +oo] comply with Assumption 13.91 and with the following chain rule: for every u G 
W^'HO.TjX) and£€ V-^T-.X*) such that 

sup £* (u(i)) < oo, G dE t (u(t)) for a.a. t G (0, T) , f f a (u(t), -£(*)) dt < oo, 

te[o,T] Jo 

then 

(du) sin = =► (cLB) s i n = (4.9) 

and the chain rule inequality (|3l£zij) holds. 

Let (u,£,E) G W([0,T];X) x L^O.TjX*) x W([0,T]) fulfill @T7]) and ([4T5]). Suppose in 
addition that u G W ' (0, T; X). Then, 

EeW hl (0,T). (4.10) 

Furthermore, the pair (u, £) fulfills 

- G a(u(t)) for a.a.t G (0,T), (4.11) 
and t/iere ZioWs £/ie improved energy inequality 



(4.12) 



Finally, if £ afeo fulfills the enhanced chain rule (|3.19l) . £/ien (|4.12l) ZioZds as an equality /or 
every < s < t < T. 

Proof. Since u G W ' (0, T; X), its distributional derivative du has zero singular part, viz. du = 
UgjcL. Then, it follows from (|4.9[) that dE = E ac L, viz. E is absolutely continuous. Therefore, 
(|4.8p becomes 

B(t) + / a («(*), -e(i)) = ft£t(«(«)) for a.a.t G (0,T). (4.13) 
Now, combining this with the chain rule inequality (!31£ztP . we conclude that f a (u(t), —£(£)) < 
(-£(t),u(i)) for almost all t G (0,T), hence (|4TT|) holds. Then, to prove (|jTT2|) we integrate 
()4.13|) . thus obtaining 

E(t) + [ f a (u(r), -£(r)) dr = £(s) + / <9 f £ r (u(r)) dr for all < s < i < T, (4.14) 



£t(u(t)) + y fa (i(r), -f(r)) dr 

< E s (u(s)) + / <9 4 £ r (n(r))dr /or aJlte (0,T], /or a.a. s G (0,i) and /or s = 0. 
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and we use 14.70 . 

If moreover £ complies with the chain rule (13.190 . then E(t) = &t(u(t)) for all t £ [0, T], since 
both functions t h-> E(t) and i H > Et(u(t)) are continuous on [0,T] and coincide on a set of full 
Lebesgue measure. Therefore from (|4.14|) we get (|4.12j) for i n> £ t (u(i)). This concludes the 
proof. □ 

As a straightforward consequence of Prop. 14.71 we have the following result, showing that, 
under a stronger coercivity assumption on the sequence of maximal monotone operators (a n ) (cf. 
(|4.15j) below), any sequence (u n ) of solutions to (|4.2j) converges up to a subsequence to a curve 
complying with (|4.10[) - (|4.12[) . In particular, observe that, unlike in (|3lai j) . in (|4.15j) we do not 
allow the "degenerate" value 1 for exponent p. Indeed, Theorem 14.81 below for instance applies 
to a sequence of operators a n — dip n , with ip n (v) — l/p n ||u|| Pn and p„ i p > 1 as n —> oo. In 
this way, we obtain a stability result for doubly nonlinear differential inclusions driven by viscous 
dissipation potentials, which generalizes the results in (TJ Thms. 3.1, 3.2]. 

Theorem 4.8. In the frame of (|3. 151) . suppose that £ : [0,T] x X — > (— oo, +oo] complies with 
Assumptions 13.91 and 14.41 Let a n : X =4 A* fulfill (l3lo;nl) /or ol/ n 6 N and 

3ci,c 2 ,c 3 >0, p > 1, <?>1 VneN V(i,!/)ea„ : 

M^llzf + c^-cs. (4 ' 15) 

Suppose that there exists a : X =4 -X - * fulfilling (|3lan[) suc/i f/iai a„ — — » a. Let (uq) <Z D be a 
sequence of initial data fulfilling (|4.3j) and Zet (w n ,£n) C W ' (0,T;X) x L (0,T; JT*) 6e solutions 
to (|4~2|) . fulfilling u n (0) = Uq and (|3.21|) /or ewery n £ N. 

Then, there exists u £ W ,p (0, T; X) tuif/i u(0) = no smc/i i/iai up to a (noi relabeled) subsequence 

u n {t) ^ u{t) for allt e[Q,T], u n ->> u inW 1<p (0,T;X), (4.16) 

and there exists £ € L 9 (0,T; A*) snc/i £/ia£ the pair (u,£) is a solution to (|3.17[) m i/ie sense of 
Definition \3.11[ fulfilling the improved energy inequality (|4. 12[) . 

The proo/ is outlined at the end of Sec. [SJ 

The rate-independent case. Let us now focus on the case in which 

a = dtp with tp : X — > [0, +oo) convex, lower semicontinuous and 1-positively homogeneous 

(4.17) 

with associated polar set A* C X* . In this case, the energy inequality (|4.5j) rephrases in a more 
explicit way. 

Proposition 4.9. Assume f3T5]l . Let a fulfill (|4"T7|l and Zei G W(0,T; X) xL^O.T : X*) 

satisfy the energy inequality ()4.5[l . Then, (u,£) fulfill 

- f (i) G A"* /or a.a. i e (0, T), (4.18) 

£ t (u(*)) + Var v ,(u; [0,t]) < £ («(0)) + / <9 t £ s (u(s))ds /or a// i £ [0, T] , (4.19) 

Jo 

with Var^ /rom (13.61) . 

In the frame of rate-independent evolution, ([4.180 is interpreted as a local stability condition, 
while the energy inequality (|4.19l) balances the stored energy £ t (u(i)) and the dissipated energy 
Var^(w; [0, t]), with the initial energy and the work of the external forces J* dt£. s (u(s)) ds. In fact, 
the local stability (|4.18[) and the energy inequality (|4.19j) yield (a slightly weaker version of) the 
notion of local solution to the rate- independent system (X,E,ip) from [49, 50 . Therein, it was 
observed that this concept is the weakest among all notions of rate-independent evolution, in that 
it yields the least precise information on the behavior of the solution at jump points. On the other 
hand, local solutions arise in the limit of a very broad class of approximations of rate-independent 
systems. This is in the same spirit as the stability results of this work. In particular, notice that 
the maximal monotone operators a n converging in the sense of graphs to a = dtp need not be 



STABILITY RESULTS FOR DOUBLY-NONLINEAR EQUATIONS 



17 



cyclically monotone. An example in this direction in the plane X = M 2 is given by the graphs 
a n — dip + (l/n)Q where Q is a rotation of tt/2. In this case a n — ^ dip but each a n is noncyclic. 



We now proceed with the 

Proof of Proposition \J^[ Let ( u,£) 6 BV(0,T;X) x L^O,! 1 : X*) fulfill (|45j) . Now, in view of 
Proposition 12.51 and of formula (|2.12b[) . we have 

/«(««(*), =^(«ac(*)) + V'*(-f(*)) = ^(«ac(t))+ijr-(-f(*)) for a.a.i e (0,T). (4.20) 
Furthermore, we have that 

/"(urinCtJ.O) > sup ($,u sin (t))=ip(u sin (t)) for ||(du) sin ||-a.a. i e (0,T) (4.21) 

where the first inequality is due to p. 141) and the second identity to (|2.12d|) . 
Then, taking into account formula (|3.9p for Var^, (|4.5p yields 

8 t («(t))+Var^(tt;[Q,i])+ / /*.(-£(«)) da < £ (u(0)) + / 9 t £ s (u(s))ds for all i S [0, T], 

Jo Jo 

which is equivalent to (|4.18p - (|4.19|) . □ 

4.2. Sufficient conditions for closedness and chain rule. Following [ST], we now show that 
conditions of A-convexity type on the energy functional £ ensure the validity of the closedness 
property (|3l£qp . of the chain rules (|3.19j) . (|3I£a|) . and of property (|4.9[l . 

More precisely, in [511 Sec. 2] the following subdifferentiability property was introduced. 

Definition 4.10. Let £ : [0,T] x X -> (-co, +oo] fulfill (|3T£^) . For every i? > 0, set 

£>« = {it £ I? : 9(u) < R} . 

We say that £ is uniformly subdifferentiable (w.r.t. the variable u) if for all R > there exists a 
modulus of subdifferentiability uj r : [0, T]x DrX Dr — > [0, +oo) such that for all t € [0, T]: 

wf^u, it) = for every it £ 
the map (t,u,v) i— > uj^(u,v) is upper semicontinuous, and (4.22) 
£t(u) — £t(u) — (£, w — u) > —u)^(u, v)\\v — u|| for all u, v € Dr and £ € <9£ t (tt). 

It was shown in [STJ Sec. 2] that, a sufficient condition for (I4.22[) is that the map u i-> £f(tt) is 
A- convex uniformly in t G [0, T], namely 

BAel V< e [0,T] Vu , ui e D V0 e [0, 1] : 

£ t ((l - 8)u + 6 Ul ) < (1 - 0)£ t (uo) + 9£ f (ui) - ^0(1 - 0)||u o - "ill 2 • 

Suitable perturbations of A-convex functionals also fulfill the closedness and the chain rule prop- 
erties: we refer to [661 1511 165] for more details and explicit examples. 

We have the following 

Proposition 4.11. Let £:[0,T|xI-) (— oo,+oo] fulfill (|31£nl) . (|3I£t|) . and the uniform subd- 
ifferentiability condition (I4.22[) . Then, £ complies with the closedness condition (|31£^p . with the 
chain rules (|3.19[) and (|51£ztP . and with property (I4.0[) . 

Proof. In (STJ Prop. 2.4], it was proved that condition P~2"2l implies (|31£^P and (|3"7T9"|) . The validity 
of ([31£^P and (|4.9I) can be checked trivially adapting the arguments developed for the proof of [STJ 
Prop. 2.4], to which the reader is referred. □ 
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4.3. Examples of quasistatic limits. Our approach to the approximation of doubly nonlinear 
evolution equations in particular allows us to discuss quasistatic limits of dynamical problems. 
Indeed, the flexibility in the choice of the approximating graphs a n , possibly noncyclic monotone, 
makes it possible to take rate-independent limits of Hamiltonian systems. We shall provide here 
some examples of ODEs and PDEs that can be reformulated within our frame. 

Let us start by considering the case of a nonlinearly damped oscillator. In particular, let q = q(t) 
represent the set of generalized coordinates of the system, M be the mass matrix, and U = U(q) 
its smooth and coercive potential energy. Assume moreover that the system dissipates energy 
in terms of a positively 1-homogeneous and nondegenerate dissipation potential D = D(q). By 
rescaling time t as et, the quasistatic limit of the system corresponds to the limit as e — > in the 
equation 

e 2 Mq + dD(q) + VU{q) 3 0. (4.24) 

The latter can be rephrased as a single doubly nonlinear Hamiltonian system in the pair v = (p, q), 
by introducing the Hamiltonian H(p,q) = U(p) + q-M q/2, the symplectic operator 




and the dissipation potential D(p,q) = D(p). Then, (|4.24|) reads 

dD(p, q) + eJ(p, q) + VH(p, q) 3 (0, 0), (4.25) 
which can be equivalently rewritten as 

dD(p) +eq + VU(p) 3 0, 
-ep + M^q = 0, 

Taking the quasistatic limit e — »• in relation (|4.25[) requires to deal with the graphs a e = dD+eJ, 
which are noncyclic monotone for all e > 0. Apart from the coercivity assumption (|3lai [) (which 
can however be relaxed in this case), this situation fits into our theory. In particular, solution 
trajectories to the dynamic problem (|4.24l) converge to solutions of the corresponding quasistatic 
limit. By generalizing the choice of the graphs a e , convergence can be obtained for a large class 
of different approximating problems. 

The nonlinear oscillator example can be turned into a first PDE example by considering the 
nonlinearly damped semilinear wave equation 

e 2 u t t + dD{u t ) - Au + f(u) = 0. (4.26) 

This is to be posed in the cylinder £1 x (0,T) for some smoothly bounded open set f2 C M. n , 
along with the positively 1-homogeneous and nondegenerate dissipation potential D, the smooth 
and polynomially bounded function /, and suitable initial and homogeneous Dirichlet boundary 
conditions (for simplicity). Equation (|4.26j) can be variationally reformulated in terms of a first- 
order system as 

&D(ut,v t )+e3(ut,v t ) + dK(u,v)3(Q,Q) inlTxV* for a.a. t e (0, T), (4.27) 

where U = H^(Cl), V = L 2 (n), the functionals D : V 2 -> [0,oo], and IK : UxV ->■ (-00,00] are 
given by 

T>(ut,v t ) = D(u t )dx, %(u,v) = J (^\Vu\ 2 + f(u) + ^\v\ 2 ^j dx 

for /' = /, and 3(ut,vt)(x) — J(ut(x),vt(x)) for almost every x G fi. Equation (|4.27j) fits in 
our frame along with the choice a e = 9D + ed, which are noncyclic monotone for all e > 0. In 
particular, owing to our analysis we can take the quasistatic limit e — > in the latter (again by 
suitably circumventing the lack of coercivity, which is inessential here) . 

Let us now provide a second PDE example by considering the quasistatic limit in linearized 
elastoplasticity with linear kinematic hardening |40) . We let Q C K 3 be the reference configura- 
tion of an elastoplastic body which is subject to a displacement u : Q, -3- M 3 and a plastic strain 
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p : f2 — »• K^ v 3 (traceless or deviatoric symmetric 3x3 tensors). Then, the evolution of the elasto- 
plastic medium is described by the system of the (time-rescaled) momentum balance (in M 3 ) and 
constitutive equation (in K^f) as 

e 2 P u tt -V{C(e(u)-p)) = b, 
dD{p t )+Wp = <C{e{u)-p) 

in D, x (0, T), where p — p(x) stands for the material density, C is the elasticity tensor (symmetric, 
positive definite), e(u) = (Vu+Vit T )/2 is the symmetrized strain gradient, b = b(t,x) denotes 
some body force density, H is the hardening tensor, and D is a positively 1-homogeneous and 
nondegenerate dissipation potential. The choice D(p t ) = R\pt\ for some R > corresponds to 
the classical Von Mises plasticity. We shall close the latter elastoplasticity system by imposing 
homogeneous Dirichlet conditions on u and no-traction conditions at the boundary (for simplicity). 
Then, the system can be recast in the form of a first-order system by augmenting the variables, 
including the momentum vt — pu t . In particular, we can variationally reformulate the system as 

dl>(u t ,v t ,p t ) + e3(ut,v t ,p t ) +d'K(u,v,p) 3 (6,0,0) mU*xV*xT for a.a. t € (0, T), (4.28) 

where now the spaces are defined as It = {u € i/^fijR 3 ) : u = in dfl}, V = L 2 (0;K 3 ), 
CP = L 2 (Q: Kjc V 3 )- The functionals and the operator are given by 

V(u t ,v t ,p t ) = f D( Pt )dx V Pt GiH^O- 
Jn 

1, , , , „, , , s 1 „„ 1 



*"(»■ ''■/') = / ( ^(e(u)-p);C(s(u)-p) + ^ P Mp+j-\v\ 2 ) dx V(u,«,p) £ UxVxT, 



m 

(v* 

3(u t ,v t ,p t ) = I -u t | V(//. r. P ) : u <VxT. 

Once again, the operators a e = dT> + s3 are noncyclic monotone for all e > 0. In particular, our 
analysis is suited in order to analyze the quasistatic limit e — > in the elastoplastic system (|4.28[) . 
This clearly distinguishes our frame from the former variational principle from [73] , which is of no 
use in the dynamical case. 



5. Proof of Theorem 14.51 

Outline. Our starting point is the fact that, the functions (u n ,£ n ) fulfill for every n £ N the 
energy identity 

£tM<))+ / /«„ (« n (s), -£ n (s)) ds = £ «)+ f d t E s (u n (s))ds foralUe (0,T]. (5.1) 
Jo Jo 

From (|5.1I) . we will deduce a priori estimates on the sequence (u n ,^ n ). Relying on well-known 
strong and weak compactness results, we will then prove the convergence (up to a subsequence) 
of (u n ,£ n ) to a limit pair (u,£). Hence we will pass to the limit as n — > oo in (|5.1[) . following 
the lines of the proof of |51l Thm. 4.4]. Namely, we will combine the finite-dimensional lower 
semicontinuity theorem [351 Theorem 5.27], with tools from infinite-dimensional Young measure 
theory (see Appendix [X] for some basic recaps), and refined selection arguments mutuated from 
the proof of [51] Thm. 4.4]. Such arguments will yield the existence of a function £ G i x (0, T; X*) 
such that the pair (it,£) fulfill the energy inequality (|4.5[) . 

Notation 5.1. Hereafter we will denote by the symbols C, C various positive constants, which 
may change from line to line, only depending on known quantities and in particular independent 
of ra G N. We will also use the place-holders 



E n (t) := e t (u n (t)), P n (t) := d t E t {u n (t)). 



(5.2) 
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Step 1 - A priori estimates and compactness: It follows from (I5.1[) and (cf. also 

estimate ([3T22]) ) that E n (t) < E n (0) + d f*E n (s)ds for all t G [0,T]. Since sup neN .E„(0) < C 
by (|4.3I) . applying the Gronwall Lemma we deduce sup tg [ T ]{i?„(t) : i S [0,T]} < C. Therefore, 
in view of assumption (|31£tP and property (|3.16p . we conclude that 

3O0VneN: sup (3(«„(i)) + \P n (t)\) < C. (5.3) 
te[o,T] 

Thanks to ([3l£i j) we then infer that 

3KmX Vn G N Vi G [0,T] : u n (t) € K. (5.4) 
Then, taking into account that fa n (u n , — £„) > a.e. in (0,T) in view of (|2.9|) . f|5. 1|> yields 

3C>0VnGN: ||/ Q „ (ii n , -£„)|| l i (0 ,t) < C. (5.5) 
In view of assumption (|3lo;i p , from (|5.5p we conclude 

T 

ci||w n (s)||+c 2 ||^( S )||*ds<C. 

Also due to (|5.4p . we ultimately deduce that 

3C> Vn G N : |K|| w([o,t];A) + ||£»(s)]U«(o,t ; x«) < C. (5.6) 

Furthermore, from the energy identity (|5.ip we immediately infer that, setting h n (t) :— E n (t) — 
/J P n (s) ds, there holds 

h„(t) - h n (s) = - f f a Ju n (r), -£„(r)) ds < V0 < s < i < T. 

Therefore we have Var(ft n ; [0, T]) = E n {0) - E n {T)) + / Q T P n (s) ds < C thanks to (JOJ) and f£3]) . 
Since (P n ) is uniformly bounded in L°°(0,T), we conclude that 

3C>0 VneN: Var(i?„; [0, T]) <C . (5.7) 



Estimates (|5.4[) . (|5.6p . (|5.7|) . and the Helly principle guarantee that there exists a subsequence 
(rife) and functions u G W([0, T]; X) and £ G W([0,T]) such that, as k -> oo, 

(«»»(*),£*(«»*(*))) -> inXxEforallte [0,T], (5.8) 

du nk = U nk ■ £ A du in M(0, T; X). (5.9) 

Exploiting Thm. 13. 4| we decompose du as 

du = (dw) ac + (du) sin = zi ac £ +ii s ||(dn) sin ||. 

Observe that, by the lower semicontinuity (|3l£nP . 

E(t) > E t (u(t)) for all t G [0, T]. (5.10) 

Further, in view of estimate (|5.5p . there exists \i G M(0,T) such that (up to not relabeled a 
subsequence) 

/»»(«*(■)»-&(•)) A M inM(0,T) (5.11) 
Moreover, by an infinite-dimensional version of the fundamental compactness theorem of Young 
measure theory (cf. Thm. [A~3l in Appendix |A"]) . we can associate with (possibly a subsequence of) 
(6ifc i Pn k ) a limiting Young measure (c t ) te ( ,T) € f (0,T;a"x E) such that, for almost all t G (0, T) 
it holds Ot(X xl) = l and er t is supported on the set of the limit points of (£„,. (t), P nk (t)) w.r.t. 
the weak topology on X* x R, viz. 

supp (<r t ) C f) {(U(t),Pn k (t))--k>j} WCak (5.12) 
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(where with g weak we denote the closure of a set B C X* x M w.r.t. the weak topology). Further- 
more, it holds 

U -» / C<M£p)=:£ inL«(0,T;X*) and (5.13) 

JX* xR 

P„ fc A /" pda t (C,p) =: P in L°°(0,T) . (5.14) 

Step 2 - Nonemptyness of admissible sets: From now on, for the sake of simplicity, we shall 
write k instead of rifc . There exists a negligible set N C (0, T) such that for every t £ (0,T)\N 
convergences (|5.8p and the support property (|5.12|) hold. Taking into account the closedness 
condition (|3l£c.p . it can be easily checked (cf. also [51] Sec. 6]), that for almost all t G (0,T) there 
holds 

(t,u(t)) e D{8E), 

£*(«(*)) = £(t), £<,(u(0)) = £7(0), (5.15) 
su PP (a t ) c {(C,j») e^* xR:CG0£t(«(t)), P<9t£i(w(t))}- 
In particular, from (|5.14j) and the third of ()5.15|) it follows that 

P(t) < d t E t (u(t)) for a. a. t G (0, T) . (5.16) 

Step 3 - liminf result for the Fitzpatrick function: In the next lines, we are going to prove 
that 

liminf / f ak (v,k(r), -&(r)) dr 

k— ¥oo / n 

t t ( 5 - 17 ) 

>[ [ fa«o(r),-0 da r (C,p)dr+ [ f™ (u s (r),0) \\(du) sin \\(r) , 

JO JX'xR JO 

In order to do so, employing |35[ Corollary 1.116, p. 75], we decompose the measure /i from (|5.11|) 
as follows: there exist /i ac , p, s i n , fj,j_ in M(0,T) such that 

A*ac < ||(dit) ac || , ^ sin < ||(du) sin || , fi± ± ||(d«)ac|| + ||(du) si „|| and , . 

A* = A«ac + Msin + M-L 

In particular, /i ac is absolutely continuous w.r.t. the Lebesgue measure L. Since f ak (uk, — £fc) > 
a.e. in (0, T), we obtain /i± > 0. We will split the proof of (|5 . 1 7[) in two steps. 

First step: Now, it follows from (|3.8[) and (|5.18[) and the Radon-Nikodym property of X that the 
set of the points to G (0,T) such that <7t (X* xK) = l and 

. (d M )([^- g ,t + g]n[o,r]) 

u a c(,ioj — lim , 

e->-0 £ 
i fto+e/2 

i(to) = lim - / £ (t) dt, and (5.19) 

£->0 £ Jt _ e j2 

— — (t ) = hm < oo 

d£ £ 

has full Lebesgue measure. From now on, we shall use the notation 

Qe(to) ■= [h - e, t + e] n [0, T] with t G (0, T) such that QSTtty holds. (5.20) 
We then prove that 

d/i 



(*o)>/ /a(i«(to),-C)cK(C,jO (5.21) 

with to £ (0,? 1 ) such that (|5.19[) holds. For any such to, it is also possible to choose a vanishing 
sequence (e m ) m such that for all m G N there holds 

A* ({to - E m , ^0 + ej n [0, T]) = ( du) ({to - s m , to + e m } n [0, T]) = . (5.22) 

In order to show (|5.2ip , we will use (|2.7p . which yields 

fa{uUto),-0 = ^M{x\uUto))-{<:,x)-p a - 1 {x\x) : (or, x*) G X x X*} . (5.23) 
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In view of (I5.23|) . we thus confine ourselves to showing that 

<l/ ' ac (i ) > (^*,Mac(<o)} - / dcr to (£,p) -p a -i(x,x*) 



Jx-xR (5.24) 
for all (x, x*) G X x X* with p a -i(x, x*) < oo. 

Now, to check (|5.24[) we observe that, since — — > a also 1 — > a -1 in the graph sense, we can 
apply Theorem 12.101 to p a -i — f* . Therefore, for any (x,x*) G X x X* there exists a sequence 
(xfc,2;£) fcgN such that 

(x k ,x* k ) -> and limsupp a -i (a£,Xfc) < p a -i(ar*,a;). (5.25) 

n— ^oo k 

Combining (|5.1ip with the third of (|5.19p (for the sequence (e m ) m fulfilling (15.22p ), we have that 



d£ 



<t ) 

lim lim e" 1 / f ah (u k (t),-£ k (t))dt 



m— >oo k- 



■ x 



Qe m (*o) 



(5.26) 



> lim inf lim inf \ ( (x%, u k (t)) + (-fk(t), x k ) - p „-i (x%, x k ) ) dt, 

ro-i-oo fc^oo JQ em {to) V fc J 



where in the latter inequality we have plugged in the sequence (x k ,x* k ) from (|5.25p and applied 
formula (|5.23[) for f ak . On account of convergences (|5.9p and (|5.13p . and of the fact that (x k ,x k ) —> 
(x, x*), we have for every m G N 



[ (x* k ,u k (t)) dt ^ [ (x*,du(t)}. (5.27) 
(£k(t),x k ) dt -> / (f(t),a:\dt 



(5.28) 
(C,z) do- t (C,p) ) dt. 



'Qe m (t ) \JX*XR 

Inserting (|5~?7) - (108"|) into (1Q6"|) and using (IOS1 . we thus get 

d/Xa 



d£ 



: (*o) 



>liminf— [/ (x*,du(t))+ [ (f (~C,x) da t (£,p) -p a -i(x,x*) ) dt 

and in view of (|5.19l) we infer (I5.24p . whence the desired (|5.21l) . 
Second step: choose to G (0, T) such that it satisfies 

{du)([t -E,t +e]n[0,T}) 



Umn(h) = lim 



o ||(du)dn|l([*o-e,*o + e]n[0 > r]) ' 
° = lim „/^L M ° + £]n[ ?' r " n , and (5.29) 

d/i si n (f , _ M([to-e,*o + g]n[o,r]) 

||(d M ) sin || w iTo||(d u ) sin ||([< -£^o + £]n[o,T]) <0 °- 

The set of all io failing any of (|5.29p is a |( du) S i n ||-null set. We are now going to prove that 

n , d /T „ (*o) > / Q °°(^i„(io),0). (5.30) 

ll( dw )sin|| 

for any to G (0, T) complying with (|5.29l) . As before, we will use the notation (|5.20l) for the 
set Q £ (to) with any such to, and we choose correspondingly a vanishing sequence (e rn ) such that 
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(I5.22j) is satisfied. In order to show (|5.30l) . in view of the representation formula (|3 . 1 3[) for /£° it 
is sufficient to show that 

ii / i^n'" ii (to) > (x* , Usinito)} for a U ( x , x *) such that p^-ifa;, a;*) < oo. (5.31) 
IKduJsinll 

With the same argument as in the previous lines, we see that 
d/x s 



|(du) sln || 



(*o) 



= I™ Km j-— - — — -— - / /„„ (u k (t), (*)) dt 

m^ocfc^oo ||(du) sill || {Q Sm {t )) jQ Bm (t ) 

> liminf liminf — — — J— — — - / ((x%,u k (t)) + (-&(t),x fc > - p ah -x(x* k> x k )) dt , 

ro->oo fe^oo ||(du) sill || (Qe m (to)) JQ Cm (t ) 

(5.32) 

where (x k ,x k ) as in (|5.25[) approximates (x,x*) from (|5.31|) . Once again, due to ()5.9[) and (|5.13|) 
we have for every fixed m £ N that 



{x* k ,ii k (t)) dt ->• / (a;*, du(t)) and 

Qs m (*o) JQe m {to) 

(tk(t),x k ) dt -> / dt. 
By construction (cf. (|5.29p ). there holds 

lim ttt^ — rTTTi — TTvT / (x*, d«(i)) = {x*,M sin (* )) 

m-i-oo ||(du) sin || (Q em (t )) JQs m (to) 
lim 7773-3 — ^ — 3-77 / (-|(t), x\ df = 0, 



m-K» ||(du) sin || (Q em (* )) JQ em (t ) 
r^oo\\(duUHQ em (t )) Pa - l[X ' X > 

We thus conclude (|5.3II) . whence (I5.30|) . 

In conclusion, passing to the limit as — >• oo in (|5 . 1 [) and relying on the initial data convergence 
(|4.3p . the energy convergence (15.81) joint with (|5.10[) . (|5.14[) . and the lower semicontinuity (I5.17p . 
we have obtained 

£*(«(*))+ T / /a Kc(s), -C) da s (Cp)ds+ [* f™(u shl (s),0)\\(du) sin \\( S ) 



JX*xI 



<£?(*) + I I U(u ac (s),-C) da s (C, P )ds+ f f™(u sin (s),0)\\(du) sin \\(s) (5.33) 

JO Jx*xM JO 



< £ (u(0)) + 



(( pda s (C, p) ds for all t G (0, T] 



Step 4 - Enhanced support properties of the Young measure (er f ) te ( ,T) : We can now 
improve the third of f)5 . 15[) . showing that indeed 

supp (a t ) C {(C,p) a*xl:(e 0£ t (u(t)), -£ G a(«(t)) , p < d 4 £ t (u(t))} for a.a.t G (0,T). 

(5.34) 

To this end, observe that, passing to the limit as n k — > oo in (|5.1j) (written on the interval (s,t)), 
yields, in view of convergences (|5.8p . (|5.1ip . and (|5.14[) . the following energy identity 

E{t) + fi([s, t}) = E(s) + J p(r) dr for all < s < t < T with (5.35a) 
H([s,t])> I I /a(6«(r),-C) da r (C,p)dr (5.35b) 
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the latter inequality due to (|5.21|) . In particular, observe that 

/«(««(*)> ~0 da t (C,p)dt < oo. (5.36) 



o Jx 



Let T C [0, T] be the set of all Lebesgue points t of P (|5.14|) . such that relations (j5.19[) and (|5.21|) 
hold, and 



E ac (t ) = lim 



S (to + |) — (to - §, 



Then T has full measure. Let us now choose a sequence (e m ), e m 4 0, such that (|5.22|) holds. 
Then, passing to the limit as e m i in (15.35)) (written for s — to — s m /2 and t = to + s m /2), we 
obtain 

E ac (t a )+ [ /a(«ac(*o).-C) do-t,, (C,p) < / pd<7t (C,p) for aU to 6 0* (5.37) 

(up to removing from T a set of zero Lebesgue measure). Now, observe that thanks to the third 
of (|5.15|) and (|5.36j) . the Young measure (crt)tE(Q,T) satisfies the assumptions of the forthcoming 
Lemma lA.41 Therefore, in view of the Young measure version of the chain rule inequality plS^j) 
therein, we find that 



-£ ac (to)+ / pdcr to ((,p)< (-C,«ac(*o)) do-to (C.P) for almost all to eT. (5.38) 

JX*xR JX'xl 

Combining (|5.37l) and (|5.38p . we deduce that for almost all to G 7 (and hence for almost all 
to € (0,T)) it holds 



(/a(«ac(*o), -0 ~ (-C,«ac(*o)» d<T to (C,i>) <0. (5.39) 

IX* xl 

Since f a is a representative function for a, we easily see that (|5.39[) holds as an equality, and that 
in fact 

-C 6 a(uac(to)) for ot -a.a. (C,p) £ supp(o to ). 

Since to £ (0,T) is arbitrary out of a Lebesgue- null set, we have ultimately proved the desired 
support property (15.341) . Furthermore, as a by-product of (|5.37[) - (15.39[) holding as equalities, we 
infer the following pointwise energy equality 



Sac(t)+/ /a(uac(t),-C)<MC,P)= / P^t(C,p) for a.a. t £ (0, T). (5.40) 

Step 5 - Selection argument and conclusion of the proof: We can now apply Lemma [A. 5 1 
and deduce that there exist measurable functions £ : (0,T) — > X and P : (0,T) — »• R such that 

(£(f),p(t)) £ argmin{/ a (u ac (t),-C) -p: (£p) £ § (t, u(t), ««(*))} for a.a. t G (0,T), (5.41) 

with §(t,u(t),u ac (t)) := {(C,p) G X* xl : ( £ 0£ t (u(t)), -£ £ a(u ac (t)), p < $£ t (u(t))}. In 
particular, 

£(t) G 0£ t (u(t)) , -£(t) G a(u«(t)), and P(t) < d t £ t (u(t)) for a.a. t G (0, T), (5.42) 
We then have the following chain of inequalities for almost all t G (0, T) 



-E ac (t) = / (/a(Uac(t),-C) -P) da t (C,p) 

Jx*xR 

>/a(«ac(t),-£(t))-P(t) ( 5 - 43 ) 
> (-f(t),U BC (i))-ft8 t (u(t)) > -Pac(t), 

where the first identity follows from (I5.40[) . the second inequality from (|5.41l) . the third one from 
the fact that f a is a representative function for a and from (|5.42[) . and the last one from the chain 
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rule inequality (|3l£zi[) . Therefore we infer that all inequalities in (|5.43[) hold as equalities, which 
proves (|4.8I) . In particular, we have that for almost all t G (0,T) 



P(t) = / pda t ((,p) = d t £ t (u(t)), 

JX'xR 

Combining (|5.44p with (|5.33[) we ultimately deduce (|4~5i Now, from (|4~5j) with (j4~6f we have that 
/ Q T | (— u ac (i)) di < oo. Since a complies with (|4.1|) . we then conclude that £ G £ 9 (0, T; X*). 
This completes the proof. □ 

Remark 5.2 (The role of the Fitzpatrick function). As pointed out in Remark l3.131 the variational 
reformulation of the doubly nonlinear differential inclusion (|3 . 1 T[) could be given in terms of any 
representative functional for a. The distinguished role of the Fitzpatrick function f a is apparent 
in the passage to the limit argument developed in Step 3 of the proof of Thm. 14.51 Therein (cf. 
(|5T2Tj) - (f5T2"g)) '), we exploit the duality formula (j2~7f for f a , as well as Theorem I2TTU1 

Remark 5.3 (Refinement of the measurable selection argument). A close perusal of Step 5 in the 
above proof reveals that, in principle, it should be sufficient to select the functions t (£(i),p(t)) in 
the set S(t, u(t)) :— {{(,p) G X* xR : ( e d£. t (u(t)), p < <9 t £ t (u(i))}, i.e. dropping the requirement 
— £ G a(u ac (t)). Indeed, if we were in the position of applying Lemma lA.51 to the set S, from the 
chain of inequalities (|5.43l) the second of (|5.44l) would still follow, yielding — G a(u ac (i)) for 
almost all t G (0,T), i.e. filfl) . 

Nonetheless, the extension of Lemma IA.5I to the set § seems to be an open problem, at the 
moment, cf. the upcoming Remark I A. 61 

We conclude this section with the 



Proof of Theorem \4-8\ Repeating the calculations from Step 1 of the proof of Thm. 14.51 we prove 
that the sequence (u n ) is bounded in W /1,p (0,T; X) and in addition fulfills estimate (I5.3p . There- 
fore, convergence (|4.16[) holds. We use the arguments from the above Steps 1 and 4 to in- 
fer that there exist G L^CLTjX*) x BV([0,T}) complying with CO} and g^). Since 
u G W 1,P (Q,T; X), Proposition 14. 71 applies, and we conclude the proof. □ 



Appendix A. Young measure results 

We fix here some definitions and results on parameterized (or Young) measures (see e.g. [T2J, f!3l 
1141177] ) with values in a reflexive Banach space y. In particular, in Section[5]the upcoming results 
are applied to the space y = X* x E. 

Notation A.l. In what follows, we will denote by ^(o,t) the cr-algebra of the Lebesgue measurable 
subsets of (0,T) and by the Borel cr-algebra of y. We use the symbol (g> for product o- 

algebrae. We recall that a Jz?(o,t) ® -measurable function h : (0, T) x y — > (— oo,+oo] is a 
normal integrand if for a. a. t G (0,T) the map y i— > h t (y) = h(t,y) is lower semicontinuous on y. 

We consider the space y endowed with the weak topology, and say that a -S?(o,t) ® ^(y)~ 
measurable functional h : (0, T) x y — > (— oo, +oo] is a weakly-normal integrand if for a. a. t G (0, T) 
the map 

y M> h(t,y) is sequentially lower semicontinuous on y w.r.t. the weak topology. (A.l) 

We denote by ^#(0,T;y) the set of all «Sf( ,T)- measura ble functions y : (0,T) — )• y. A sequence 
(y n ) C ^#(0,T;y) is said to be weakly-tight if there exists a weakly-normal integrand /i : (0,T) x 
y — > (— oo, +oo] such that the map 

y i— > h t (y) has compact sublevels w.r.t. the weak topology of y , and sup / h(t, y n (t j) dt < oo. 

n Jo 
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Definition A. 2 (Young measures with values in y). A (time-dependent) Young measure 
in the space y is a family cr := {&t}te(o.T) of Borel probability measures on y parameterized by 
t G (0, T), such that the map on (0, T) 

th^at(B) is Jf( ,T)-measurable for all B G ®{$). (A.2) 

We denote by W (0, T; y) the set of all Young measures in y. 

The following result is taken from [51] (cf. Thms. A.2 and A.3 therein). It is a generalization of 
the so-called Fundamental Theorem of Young measures (cf. the classical results |121 Thm. 1], [131 
Thm. 2.2], [14], [77] Thm. 16]), to the case of Young measures with values in y endowed with the 
weak topology (see also [661 Thm. 3.2] for the case in which y is a Hilbert space endowed with the 
weak topology). 

Theorem A.3 (The Fundamental Theorem for weak topologies). Let (y n ) C ^#(0,T;y) be a 
weakly-tight sequence. Then, 

(1) there exists a subsequence (y nk ) and a Young measure cr = (ct) tg / T) ^ ^(0, T;y) such 
that 

oo 

limsup \\y nk (t)h < 00 and supp(er t ) C Q {y nk (t) : k > j} for a.a.t G (0,T), (A.3) 

fcfoo . =1 

(where B denotes the closure of a set B C y w.r.t. the weak topology), and such that 
for every weakly-normal integrand h : [0, T] x y — > (— oo,oo] suc/i i/ia£ /i _ (•, y„ fc (•)) is 
uniformly integrable it holds 

cT pT 



liminf / h(t,y„ k (t)) dt> h (t,y) da t (y) dt . (A.4) 

Jo Jo Jy 

(2) In particular, let (y n ) C L q (0,T;y) be a bounded sequence, with q G (1, +oo]. Then, there 
exists a further (not relabeled) subsequence (y nk ) and a Young measure cr = {&t}te(o,T) € 
T;V) smc/i that for a. a. t G (0,T) properties (|A.3j) /jo/gL Setting y(t) := J„ y dtrt(j/) 
/or almost all t G (0,T), i/tere /to/ds 

y„ fc -y mL p (0,T;y). (A.5) 

A. l. A Young- measure version of the chain rule. In what follows, we will work with Young 
measures with values in the space y = X* x M. Our first result, a small variation of [5T] Prop. 

B. l], provides the version of the chain rule inequality (13l£zi|) in terms of Young measures used in 
Step 4 of the proof of Thm. PI 



Lemma A.4. In the frame of (|3.15[) . let a : X X* fulfill (|31anP and the coercivity condition 
P~T]) . and let £ : [0,T] x X -> (-oo, +oo] comply with Assumvtion WM Let u G BV([0,T];X) 
satisfy 

sup £ t (u(t)) < oo, (t,u[t)) G dom(d£) for a.a. t G (0,T), / |d t £ t (it(£))| di < oo, 
te[o,T] Jo (A.6) 

G W([0, T]) smc/i that E{t) = E t (u(t)) for a.a. t G (0, T), 
and let (o~t)te{a,T) G £9^(0, T; A* xl) 6e a Young measure such that 

V(f,p) G supp(a t ) : £ G 0£*(u(t)) , p < $£*(«(«)) for a.a. t G (0,T) (A.7) 

/a(Uac(s), -C) d(T(C,p)ds < OO (A.- 



Then, for almost all t G (0, T) suc/i that t is Lebesgue point of E ac and u ac there holds 

E ac (t)> [ ((C,u ac (t))+p)da t (C,p). (A.9) 

JX'xR 
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Proof. We consider the set K(t,u(t)) := {(£,p) G A* x R : £ G dE. t (u(t)), p < d t £. t (u(t))} . 
Repeating the very same arguments as in the proof of [511 Prop. B.l], we can show that there 
exists a sequence (£ n ,Pn) of strongly measurable functions (£, n ,Pn) ■ (0,T) — » X* x R such that 



{(&,(*),?»(*)) : n G N} C K(t,u(t)) C {(f n (*),Pn(*)) : n G N } for a. a. t G (0,T) (A.10) 
(where £? denotes the closure of B C X* x R w.r.t. the strong topology of A* x R). 
We now claim that the sequence (^„,p„) can be chosen such that 

Vn £ N : £ n G L\0,T-X*) and sup f f a (u ac {t), -£„(<)) dt < oo . (A.ll) 

nSN Jo 

To this aim, we define the function g(t) := inf{/ Q (u ac (t), — f) '■ (CiP) G A(i,u(i))} for almost all 
t G (0,T). Notice that due to jlTOjl it holds 

g(t) := inf{/a(«»c(<), -&.(*))} for a - a - * e (°> T ) ( A - 12 ) 

and hence g is measurable on (0, T). Moreover, 

S(t)d*< / / /a(«ac(*)»-C)do-t(C,J>)dt<c». (A.13) 

JO JXxR 

With a straightforward adaptation of the argument of [5TJ Prop. B.l] (see also [551 Lemma 3.4]), 
from (|A~12l) and (lA~13j) we deduce (|A~TT|) . 

In view of the obtained (|A.10[) and (lA.lip . we are in the position to apply the chain rule 
inequality (|31£zi[) to the pair («,£„) for every n G N. Therefore for every n G N there exists a set 
T„ C (0,T) of full measure such that E ac {t) > (£ n (i), u ac (t)) + p n {t) for all t G T„, where we have 
also used that p n (i) < dt&t(u(t)). The set T = r)neN^»> nas s ^ n full measure, and there holds 

for all t G 7 

Sac(*)>(C.«ac(*))+P for all (C,p) Gconv^,^)), (A.14) 
the latter set denoting the closed convex hull of K(t,u(t)). Integrating (|A.14[) w.r.t. the measure 
tr t we obtain (|A.9[) . □ 

We conclude with the measurable selection result exploited in Step 5 of the proof of Thm. 14.51 

Lemma A. 5. In the framework of (|3.15p . let a : X =4 X* fulfill (|3la:n[) and the coercivity 
condition (|4.1[) , and Zet £ : [0,T] x A — > (— oo, +oo] comply with Assumptions 13.91 and 14.41 
Furthermore, let u G BV([0,T]; A) fulfill (|A.6p . Suppose that for almost all t G (0, T) 

S(i, «(*),««(*)) := {(C,P) G A* x R : C G 9£ t (u(f)) , ~C G a(u ac (t)) ,p < d t £ t {u{t))} ? . 

(A.15) 

Then, there exist measurable functions £ : (0, T) — > A* and P : (0,T) — > R suc/i £/ia£ 
(£(t),P(t)) eargmin -C)-p:(C,p) G S(t,u(t), /or a.o. t G (0, T) . (A.16) 

Proof. The argument follows the very same lines of [511 Lemma B.2]. First of all, we observe that 

argmin{/ Q (u ac (<),-C) -p : (C,p) G S(t, «(t), ««(<))} ^ for a.a. i G (0,T). (A.17) 

To this aim, let (Cn,Pn) C &(t,u(t),u ac (t)) be an infimizing sequence: then there exist constants 
C, C > such that for every n G N 

C > / Q (u ac (t),-C„)-P„ = {-(n,U ac {t)} - Pn > Cl ||u ac (f)||f + C 2 ]|Cn||^-C 3 -C" (A.18) 

where we have used that — C„ G u ac (t), the coercivity property (|4.ip of a, and that p„ < 
dt£t(u(t)) < C due to the fact that sup tg [ ,T] &t(u(t)) < oo and to (|3I£?[) . Therefore, we in- 
fer that sup„ gN (||Cn||* + \Pn\) < oo. Hence, there exist (C,p) G A* x R such that, up to a not 
relabeled subsequence, — 1 ( in A* and p n — > p. Thanks to the closedness condition (|3i£q[) and to 
the weak closedness of a(u ac (i)), we have (C,p) G B(t, u(t), u ac (i)). Using that £ H> f a (u ac (t), — C) 
is (sequentially) weakly-lower semicontinuous, we conclude that 

liminf (/ a (M ac (i),-C„) -p n ) > / a («ac(*)j-C) -p 
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and (|A.17|) ensues. 

Once obtained (|A.17[) , the argument for (|A.16[) is a straightforward adaptation of the proof of 
[5T1 Lemma B.2], to which we refer for all details. Let us only mention here that the existence 
of (£, P) as in (|A.16|) is a consequence of the measurable selection results [211 Cor. III. 3, Thm. 
III.6]. □ 

Remark A. 6. Let us stress that the requirement ( S a(uac(i)) in the definition (|A.15[) of the set 
S(i, u(t), u ac (t)) has a crucial role in proving that 

argmin{/ Q (it ac (t), -£) -p : ((,p) e S(t, u(t), «ac(t))} 

is nonempty. In fact, it ensures the estimates in (|A.18[) for any infimizing sequence (C, n ,p„). 
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